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Abstract 

The reheating process for the inflationary scenario is investigated 
phenomenologically. The decay of the oscillating massive inflaton field 
into light bosons is modeled after an out of equilibrium mixture of in- 
teracting fluids within the framework of irreversible thermodynamics. 
Self-consistent, analytic results for the evolution of the main macro- 
scopic magnitudes like temperature and particle number densities are 
obtained. The models for linear and quadratic decay rates are investi- 
gated in the quasiperfect regime. The linear model is shown to reheat 
very slowly while the quadratic one is shown to yield explosive parti- 
cle and entropy production. The maximum reheating temperature is 
reached much faster and its magnitude is comparable with the inflaton 
mass. 
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1 Introduction 



During an inflationary stage driven by the potential energy of a classical 
scalar field, the temperature of the universe was redshifted to almost zero. 
Thus a mechanism to raise the temperature of the universe at the end of the 
inflationary stage is required to match this scenario with the standard hot 
Big Bang cosmology 0. 

Currently it is assumed that soon after the end of inflation, the scalar field 
began to oscillate near the minimum of its effective potential. Decay of this 
scalar field through particle production and thermalization of its energy with 
a huge increase in temperature and large entropy production, occurred during 
the reheating period. It started with violent particle production through 
parametric resonance, also called the preheating process, and lead to a highly 
nonequilibrium distribution of the produced particles, subsequently relaxing 
to an equilibrium state [P, ^ H. 

In this paper we will consider the decay of the massive homogeneous 
inflaton scalar field coupled with light bosons Xy it arises in chaotic 
inflationary models. A detailed quantum mechanical analysis is quite com- 
plicated and involves taking into account nonlinear processes, backreaction 
effects, rescattering, etc 0. Further, since particles created during the pre- 
heating stage are far from equilibrium, their thermalization is achieved via 
collisional relaxation. This process of thermalization has not been under- 
stood yet from first principles 0. Here, we follow an alternative description 
of the physics, and instead of both fields and their interaction we consider 
an effective, phenomenological model in terms of an out-of-equilibrium mix- 
ture of two reacting fluids within the framework of relativistic irreversible 
thermodynamics [|, 0. The difference in the equations of state of the fluid 
components and the simultaneous deviation from the detailed balance in the 
rate equations for the interfluid reactions give rise to an entropy production 
effect that manifests itself as a reactive bulk pressure This two-fluid 

model can be described in the effective one-fluid picture, based on a single 
temperature that we identify with the temperature of the reheating process. 
One of the advantages of this approach is that it is possible to treat selfcon- 
sistently the dynamics of geometry and matter |Tl[ . 



In this phenomenological approach we ignore microscopic details and con- 
centrate on most relevant macroscopic magnitudes like particle number den- 
sity, temperature, entropy production. The net effect of the interaction term 
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between both fields is collected in the particle number rate of change, and 
the finite temperature takes into account the backreaction effect of the pro- 



duced particles. A phenomenological approach was initially proposed in [|n 
The methods used in that paper for the solution of the evolution equations 
lead to a vanishing bulk viscous contribution to the speed of sound. However 
their model depends strongly on this contribution so that their procedure of 
calculation needs to be improved. In this direction we introduce a physically 
founded perturbative scheme that yields more satisfactory results. 

The plan of this paper is as follows. Section 2 introduces the two-fluid 
magnitudes and section 3 presents main equations of the model in the one- 
temperature picture. In section 4 we introduce an expansion in powers of 
a parameter that measures de departure from perfect fluid behavior. This 
allows us to solve the system of equations order by order. In section 5 the 
model with a linear ansatz for the decay rate is solved, and we show that it 
reaches the reheating temperature only for very large times. We propose in 
section 6 a quadratic ansatz for the decay rate that is better motivated both 
from microphysical and geometrical point of view. We investigate thoroughly 
the evolution of the physical magnitudes and we find that the behavior of 
this model is very interesting. Finally, conclusions are stated in section 7. 



2 Two-fluid picture 

For the two-fluid model we assume that the energy-momentum tensor T**^ 
splits into two perfect fluid parts, 

rpik rpik _j_ rj-lik ^-j^^ 

with {A = 1, 2) 

Tf = p^u'u'' +p^h''' . (2) 

is the energy density and is the equilibrium pressure of species A. For 
simplicity we assume that both components share the same 4-velocity m*. 
The quantity h^^ is the projection tensor /i**' = g^^ + u^u^. The particle flow 
vector A^^ of species A is defined as 

N: = ny , (3) 
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where is the particle number density. We are interested in situations where 
neither the particle numbers nor the energy-momenta of the components are 
separately conserved. The balance laws for the particle numbers are 

= < + = n,T, , (4) 

where = u*j is the fluid expansion and is the rate of change of the 
number of particles of species A. There is particle production for > 
and particle decay for F^ < 0, respectively. For F^ = 0, we have separate 
particle number conservation. 

Interactions between the fluid components amount to the mutual ex- 
change of energy and momentum. Consequently, there will be no local 
energy-momentum conservation for the subsystems separately. Only the 
energy-momentum tensor of the system as a whole is conserved. 

Denoting the loss- and source-terms in the separate balances by t^, we 
write 

T% = -t\ , (5) 

implying 

p, + e{p,+p,) = ux, , (6) 

and 

{p^+p^)u-+p^^,h'^^ = -h1t\ . (7) 

All the considerations to follow will be independent of the specific structure 
of the t\ . 

Each component is governed by a separate Gibbs equation 

^.d., = d^+p,d-l, (8) 

"'A ""A 

where is the entropy per particle of species A. Using eqs.(||) and one 
finds for the time behaviour of the entropy per particle 

^aTaSa = - {p^ +pjT^ . (9) 

With nonvanishing source terms in the balances for and p^, the change 
in the entropy per particle is different from zero in general. 
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The equations of state of the fluid components are assumed to have the 
general form 

Pa=Pa{t^a^T^) (10) 

and 

Pa = Pa K^^a) , (11) 

i.e., particle number densities and temperatures are regarded as the 
basic thermodynamical variables. The temperatures of the fluids are different 
in general. 

Differentiating relation ([TT|), using the balances (H) and (|^) as well as the 
general relation 

dp A ^ Pa+Pa dp^ 
that follows from the requirement that the entropy is a state function, we 



find the following expression for the temperature behaviour [O, O, 14 



T ((^ V \ dPA/9T^ , Ua^A-'^AipA+PA) .-.^N 



The entropy flow vector 5"° is defined by 



S: = n^sX , (14) 

and the contribution of component A to the entropy production density be- 
comes 



where relation (^) has been used. 

According to eq.(^) the condition of energy- momentum conservation for 
the system as a whole, 
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(r'= + rf)^ = o, (16) 

implies 

t"^ = -tl . (17) 

There is no corresponding condition, however, for the particle number bal- 
ance as a whole. Defining the integral particle number density n as 

n = rij + 71-2 , (18) 

we have 

h + Qn = nT , (19) 

where 

nV = n,r, +n^T^ . (20) 

r is the rate by which the total particle number n changes. 
The entropy per particle is |T^ 

„ _Pa+Pa f^A 



^aTa 



(21) 



where /x^ is the chemical potential of species A. Introducing the last expres- 
sion into eq. (|l5|) yields 



Sl.,a = -f^n^r, + ^ . (22) 



rji A 
^ A 

For the total entropy production density 

S-a = S:,, + 5^% (23) 



we obtain 



S- = -^^nV -(^-^\ n,T, + {---] uX • (24) 



' 2 , 



We assume that the entropy per particle of each of the components is 
preserved. The particles decay or are produced with a fixed entropy s^. 
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This 'isentropy' condition amounts to the assumption that the particles at 
any stage are amenable to a perfect fluid description. When = 0, we get 
from (^) 

M: = (p,+Pjr, , (25) 

and combining eqs.(|l7D and (^), one has 



Uaf", = {p, +Pjr, = -Uatl = - (p, + r, , (26) 

which provides us with a relation between the rates and : 
Inserting the last relation into equation (^O]) yields 



The quantities = (p^ + Pa) /"^a the enthalpies per particle. 

With the relations (^) and (^) the entropy production density (p^) 
becomes 



Pi +Pl P2+ P2 



(29) 



We emphasize that according to the equations of state (|T^) and (|TTD the 
quantities pi, p^ and depend on T^, while Pj, p^ and depend on T^. In 
general, we have 7^ T'a- 

To complete the two-fluid model we give the thermodynamical properties 
of both fluids. If we ignore effects associated with particle creation, after 
inflation the field oscillates near the point = at a frequency given by 
its mass m. Its oscillation amplitude $ falls off as t^^ and its mean energy 
density p^ = m^$^/2 decreases as a~^, in the same way as pressureless dust 
Hence we represent the massive inflaton field as a fluid, named 1, 
of nonrelativistic particles with mass m, energy density pi = p^, pressure 
Pi <^ pi and number density rii ~ pi/m. On the other hand, the decay 
products of the scalar field are ultrarelativistic for m ^ m^, and we model 
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them as an ideal relativistic fluid of massless particles, with energy density p2, 
equilibrium pressure p2 = P2/3 and particle density n2. We also assume that 
both fluids are classical and ideal and that each species remain in thermal 
equilibrium along the reheating process. That is, we neglect dissipative effects 
for each of the fluids in comparison with particle production effects. 



3 Effective one-temperature picture 

We have also the effective one-temperature alternative description of a two- 
component fluid that is based on a single Gibbs-equation for the system as 
a whole: 

Tds = d^ + pd- - (/i, - /i J d^ , (30) 

n n n 

where p is the equilibrium pressure, p is the energy density and s is the 
entropy per particle. The temperature T is the equilibrium temperature of 
the whole system and is defined by |IB, 



pA^^^T^) + p^{n^,T^) = p{n,n^,T) . (31) 

Furthermore, we assume that the cosmic fluid as a whole is characterized by 
the equations of state 

p = p{n,n^,T) (32) 

and 



p = p{n,n^,T) , (33) 

The temperatures and T,^ do not appear as variables in the effective 
one-temperature description. An (approximate) equilibrium for the entire 
system is assumed to be established through the interactions between the 
subsystems on the right-hand side of the balances (^) and (0) . In this picture 
the cosmic fluid splits into two effective fluid components whose equations 
of state are assumed to have the same form as those of the "real" fluids, 
but shear this common temperature. So, these two effective fluids are in a 
(formal) thermal equilibrium between them while the "real" fluids are not. 
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If we set = = T in (|21|) and make the sphtting for the energy 
density of the system 



p (T) = (T) + (T) (34) 
and for its equihbrium pressure 

p(T)=p,(T)+p,(T), (35) 

the description based on relation (|30|) is consistent with the description re- 
lying on the relations (|]) for ns (T) = n^s^ (T) + n^s.^ (T). As long as the 
pressures are those for classical gases, = n^T, the equilibrium pressure 
p of the system as a whole depends on n = n^^ + only and the separate 
dependence on on the right-hand side of eq . (p2D may be omitted. 

For components out of equilibrium, even when F = = = 0, the 
total pressure is generally different to the sum of the partial pressures. In 
this case we have |T^ 

Pi ('^n^i) +P2 (^2^^^) =p{n,T) + TTd . (36) 
instead of (|53|). Here tt^^ is the viscous pressure arising from differential 



temperature variation rate between both fluids fl^. From eq.(|T3D the cooling 
rate T^/T^ is different from T^/T^ even for F^ = F2 = if the subsystems 
are governed by different equations of state. The expansion of the Universe 
tends to increase the difference between and T^. 

On the other hand, deviations from detailed balance, i.e., F^ 7^ 0, leading 
to F 7^ in general, generate an effective 'reactive' bulk pressure vr. For the 
corresponding energy-momentum tensor of the system as a whole we write 

T'' = pu'u^ + (p + tt) h'^ . (37) 

The reactive bulk pressure is determined by the consistency of the expression 
for the entropy production density obtained in the two-temperature picture 
(p9D and the expression obtained in the one-temperature picture 



= + ns (38) 

The last one arises from the entropy flow vector 5''^ = snu^ Q , where s is the 
entropy per particle and is the four-velocity of the fluid. It was estimated 
that TTd is one order of magnitude smaller than tt ||10|| . 
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From the Gibbs-equation (^) one finds for tlie cliange in tlie entropy per 
particle 



ns = -P-^T-"-^ {tl^) (r, - r,) , (39) 

Even for Sj = = we liave s 7^ in general. However, during the 
preheating stage, particle production dominates over per particle entropy 
change as the main entropy production source. 

The equilibrium temperature of the system T is the key magnitude of 
our model. On the one hand we identified it with the reheating tempera- 
ture, and on the other hand its definition ( |3TD provides a mean to account 
for backreaction effects of the created relativistic particles on their massive 
counterparts. Its evolution law is given by |]TU 



? = (r-e)^-^^ (40) 

T ^ ' dp/dT Tdp/dT ^ ^ 

The Einstein field equations for a spatially-fiat Robert son- Walker space- 
time are 

3H^ = np , (41) 

H = ~ip + P + n) , (42) 

where k is Einstein's gravitational constant and H = a/a, a is the cosmic 
scale factor. We use units such that c = 1, ks = 1 and h = 1, then k, = 
Sir/Mp, where Mp is the Planck mass. Finally a dot denotes derivative with 
respect to comoving time. Conservation of (|37| ) leads to 



p + 3if (7p + 7r) = (43) 

where 

Mt) . 1 . II (44) 

is the time dependent polytropic index. Using ( ^T]) and ( |4^ ) we get 

KIT = -3^H'^ - 2H (45) 
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When m ^ T, the equations of state for the nonrelativistic fluid 1 are: 

Pi = riim + \niT , pi = riiT (46) 
while for the relativistic fluid 2 are: 

P2 = 3n2T , p2 = n2T (47) 

For these fluids we have 

7 = 1 + ^ (48) 

riim + (^|ni + 3n2j T 

To calculate the evolution of the temperature we need also T. Using (p^) 
and the equations of state (|46|) and (0), we obtain 



where the decay rate Q = |ri| is an input to the model that must be chosen 
from a fundamental microphysical theory or from phenomenological consid- 
erations. 



4 The quasiperfect expansion 

We assume that the viscous effects, are small. If r is the mean interaction 
time of the particles of the fluid, we have that u = {tH)~^ is the number of 
interactions in an expansion time. Perfect fluid behavior occurs in the limit 
V ^ oo. Small departures from this behavior occur for large and a con- 
sistent hydrodynamical description of the fluids requires u > 1. Thus we are 
lead to assume that tH is small and we propose the following "quasiperfect" 
expansion in powers of i>~^ 

if = //o (l + ^ + • ■ ■) (50) 

vr = ^ + ^ + ... (51) 
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where tt^, will be fixed by the thermodynamical theory adopted. Inserting 
)0D and (|5l| ) in (^5]), we find the equations that determine the coefficients 



of the expansion up to first order in u ^ 

Ha + {il2)lHl = (52) 

k ^ (53) 

where we have assumed that f ^ tH. Solving these equations we get 

and 

In the particular case we choose the truncated transport equation of 
Causal Irreversible Thermodynamics for the bulk viscosity pressure 

vr + TTT = -3Ci^. (56) 

where we identify r with the relaxation time, the bulk viscosity coefficient C, 
is given by |19| 



- = c^(P + P), (57) 
r 

Cf, is the bulk viscous contribution to the speed of sound v, v"^ = + < 1, 
and Cs is the adiabatic sound speed, we find 

/.vr ^ -IchH', (58) 

Also it can be easily seen that Hq is an exact solution of equations (^TJ), p^), 
(HD, (ID and (13) provided = 0. 

To get an estimation of the physical parameters in the quasiperfect regime 
it is enough to keep calculations at zero order in and we can neglect the 



viscous terms in (^31), (EOl) and (08|). In this order of approximation the 



results are independent of the form of the transport equation. Nevertheless 
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this approach will allow us to give a reasonable description of the reheating 



process. Then (|3]) becomes 

p + 3i/o7P ~ (59) 



whose solution is 



Pit) ^ ^2 (60) 



Udtlit)) 

where po is a positive integration constant. Also (EDI) becomes 



Then, when the decay rate Q is large enough, it may overcompensate the 
adiabatic term during the initial 'preheating' stage and make the temperature 
rise. 



5 Linear reheating 



Following [11] we assume first a linear ansatz Q = (3H, with an adimensional 
constant /? > 0. 

To solve the coupled system of equations (|^), (|), ( |6lD and ( ]48| ) for the 
evolution for n, rii, n2, T and 7, we use an iterative scheme that starts from 
a fully nonrelativistic stage: rii = n, 7 = 1. In this regime 



Then (0) reduces to 



whose solution is 



T 
T 



2 
3t 




(62) 



(63) 



T(t) = Ti 



1 - 



(64) 
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where to is an arbitrary integration constant with dimension of time, and 
Ti = m/3/12. In order to describe the reheating scenario we need to choose 
to > 0. As in the previous inflationary stage the Universe "supercools" [Q], 
we assume that the inflation period ends at t ~ to when T ^ 0, and we 
take solution (|64D for t > to- Then the temperature grows monotonically 
and approaches asymptotically its maximum value Ti, which we could call 
the reheating temperature of the linear model, for t —>■ oo. We see from (^) 
that the temperature rises quite slowly. It takes a time over 30to for the 
temperature to approach at 1% of Ti. For this reason, we will not pursue 
further the consequences of the linear ansatz and in a following section we 
will propose an improved ansatz for Q which leads a much faster rise in 
temperature. 



6 Quadratic reheating 

In the previous section, we have seen that a linear ansatz for Q leads to a 
very slow process of reheating. A microscopic interaction term like gcf^x^i 
suggests an effective quadratic decay rate F ~ <l>^ ~ t~^. Here we look for 
an improved ansatz that takes into account the scalar nature of F. In this 
direction we propose Q = (3R, where R = Q [h + 2H^^ is the curvature 

scalar and (3 is a positive constant with dimension of time. 
At zero order in z/~^ we obtain 

Q = 6(3(2-^) Hi = QnHl (65) 



2 

This expression shows that Q > for 1 < 7 < 4/3 and vanishes in a radiation 
dominated era, when $ is very small. Inserting (|65|) in (|49| ) and using (|6l|) 
we obtain at zero order in 

and we solve it using the iterative scheme. 
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6.1 Nonrelativistic regime 

In this regime when the fluid is dominated by massive particles, equation 



(|66D becomes 

T + 2HoT ~ mriHl (67) 



whose general solution is 



^, , 4mri 



t 



(68) 



As in the linear model, the reheating scenario corresponds to to > 0, and we 
take the solution ( |5BD for t > t^. This solution also starts from T = 0, but 
this time it rises violently to a maximum temperature of reheating 

in a period tr—t^ = 37to/27, where tr is the time when maximum temperature 
occurs. Assuming that inflation ends at to = S/m, where 5 is a numeric 
constant, this rising period is 1.376/m. This is a time of the same order as the 
period of oscillation of the inflaton field 271 /m. After that, the temperature 
begins to fall. Though these calculations are carried out in the approximation 
T -C m, and this implies that ( ^ ) is strictly valid only for xi -C to, it also 
suggests that Tr may be of order m when ri is large enough. 

Inserting (|62D and (|65D in (^) and solving the resulting equation for the 
number density of massive particles we find 

g8ri/(3t) 

ni(t)~nio^^ (70) 

where uiq is a positive integration constant, so that ni{t) is a decreasing 
function. 

Using (|19D, (^91), (|68|) and ([70|) , we find the total number particle density 
in the first departure from the nonrelativistic initial state 



n{t) ~ — 



1 H / ax- 



(71) 

where x = t/to, and uq is a positive constant. Assuming that the reheating 
temperature is much smaller than m we obtain the approximated expression 
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n{t) ~ — 



(72) 



2no 

The expression ([70|) for the density of nonrelativistic particles is a de- 
creasing function, while the expression (^) for the total number density has 
a peak. The consistency requirement ni < n is satisfied after the time tp 
when ni(tp) = n(tp) = Up. This time is given by 



Xp 



1 + exp 



2 (a 



3a 



(73) 



where a = nio/rio, and marks in our model the starting point of the preheat- 
ing stage. We demand the temperature 



ntp)^4(|)'(xf -1)t. (74) 



3 

to be very low and this occurs when Xp is close to 1, so that a <^ 1. This 
condition corresponds to an explosive production of light bosons. In effect, 
from ( |1|), ( |7|) and ( |7g ) we get 

and just after tp this particle density grows linearly with an exponentially 
large slope 

n2{t) ^ —2 {x - Xp), x-Xp'^l (76) 

reaching a large peak density Up/a in a time Xm ~ 1 + 3a/4. After that, 
both n and n2 begin to fall rapidly as the effect of expansion dominates over 
particle production. 

In the regime of large particle production we can neglect in (^) the 
change in s, so that the entropy production density is approximately given 
by 



2t2 (t - t'jH^^ 
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and we obtain the following lower bound for the entropy density 



Thus we find that the generated entropy per massive particle S/n{tp) > s/a 
is very high. 

6.2 Intermediate regime 

Here pi ~ p2, i-e. riim ~ 3n2T so that n2 3> ni and 7 |. At the 
microscopic level this corresponds to the regime when the backreaction of 
the produced particles ceases to be negligible. We can estimate the time 
ti when this regime begins from the equation p{ti) = 2pi(tj), using (|i6|), 
(|60D , (|6^) and (ffOl). Therefore ti becomes a function of = mTi/6 and 
the intermediate regime is reached provided K > 3 In 2/8. For lower values 
of K particle production rate is not large enough to compensate the faster 
decrease of relativistic matter energy density due to cosmic expansion. On 
the other hand, K cannot be much larger than this lower bound if <^ m. 
For instance, we get Xi ~ 41 and Tr/m ^ 0.038 for K = 0.27 . We show in 
Fig. 1 the evolution of pi and p2 for this value of K, from fully nonrelativistic 
to intermediate regime. 

In this intermediate regime holds 

Ho{t) ^ ^ (79) 

-24n/(49t) 

m^4^'—^^ (80) 

96Ti/(49t) 

^12/7 (81) 

-72ri/(49t) 

nit)^nl;^'--^ (82) 
where Tq\ h^q and are positive integration constants fixed by matching 



at ti expressions (q^); dUl) (ID with their nonrelativistic counterparts 



, dTPP and (|7Tp. After the initial stage of large ultrarelativistic particle 
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creation in the nonrelativistic regime the nonrelativistic particle decay pro- 
cess slows down. All particle densities and the temperature decrease, though 
the dilution and cooling rates are smaller because of the slower cosmic expan- 
sion. The whole evolution of the temperature and particle number densities 
along the nonrelativistic and intermediate regimes is plotted in Figs. 2-4 for 
K = 0.27. 

6.3 Ultra-relativistic regime 

In this last stage of reheating, the energy density becomes dominated by the 
radiation fluid, i.e. ^ nim, 7 ~ |- Therefore we have 

i?. « i (83) 

rpiu) 

m ^ (84) 

- n^{t) ^ ^ (85) 

where Tq""*, Uiq and no""* are positive integration constants. This stage marks 
the end of reheating and the cosmic medium approaches a perfect, relativistic 
fluid with vanishing viscosity and particle production. The temperature and 
particle densities continue falling at about the adiabatic rate. This is smaller 
than the intermediate stage rate as the radiation-dominated universe expands 
even slower. A small remnant of nonrelativistic particles may remain. 



7 Conclusions 

With the aim of understanding the process of reheating we have developed a 
phcnomenological model of two interacting fluids. In this model the oscillat- 
ing infiaton field is modeled after a nonrelativistic fluid of massive particles 
that decay into an ultrarclativistic fluid of massless particles. 

We have carried out a "quasiperfect" perturbative expansion to solve the 
Einstein equations. This corresponds to a small viscous pressure regime. To 
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the lowest order in this expansion, the thermodynamical magnitudes do not 
depend on the transport equation for the bulk viscosity. 

A particle decay rate, linear in the expansion rate, does not yield a suit- 
able model as the maximum reheating temperature is reached only for very 
long times. However an alternative quadratic ansatz for the decay rate gives 
a good behavior of the thermodynamical magnitudes. The last one was in- 
vestigated in detail along its evolution from the initial nonrelativistic stage 
to the final domination by ultrarelativistic particles. The equilibrium tem- 
perature of the whole system and the produced particle number density rise 
violently in a time comparable to the oscillation period of the inflaton field. 
Also, large amounts of entropy are produced in this initial stage of reheating. 
The time when the energy density of the decay products becomes comparable 
with the nonrelativistic particles, is similar to the microscopic calculations 
for the time when backreaction effects become important. 

In a future paper we will investigate phenomenological reheating models 
far from the perfect fluid regime. 
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Figure Captions 



Figure 1. Plot of the equilibrium energy densities of nonrelativistic and 

relativistic particles, pi and p2 in units of po = p{to) vs. adimensional time 
X, from X = 1 to X Xi, for K — 0.27. In this case Xi ~ 41. 

Figure 2. Plot of the equilibrium temperature T in units of inflaton mass m 
vs. adimensional time x along the nonrelativistic and intermediate regimes, 
for K = 0.27. 

Figure 3. Plot of the nonrelativistic particle number density rii in units of 
the initial density vs. adimensional time x along the nonrelativistic and 
intermediate regimes, for K — 0.27. 

Figure 4. Plot of the relativistic particle number density n2 in units of 
the initial density Up vs. adimensional time x along the nonrelativistic and 
intermediate regimes, for a = 0.1 and K = 0.27. The inserted plot shows a 
detail of the evolution of n2 about its peak. 
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